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Report of the Office of 
Technology Assessment of the 
Congress of the United States 

The following text is taken from the report of the Office 
of Technology Assessment of the Congress of the 
United States, titled The Fusion Energy Program: The 
Role of TPX and Alternate Concepts. 

The Web version was dated January 31, 1995. The full 
report is available from the World-Wide Web at 
http://otabbs.ots.gov:80/E574T101. 

I have exerpted parts of Chapter 4: Alternate Concepts 
for Fusion Energy Administration. 

STATUS AND PROSPECTS OF ALTERNATE 
CONCEPTS 

There are several alternate fusion concepts with a wide- 
range of maturity levels or development of the informa- 
tion base. Over the past decades, the primary focus of 
the fusion energy program has been on several MFE 
concepts, (see footnote 9) Extensive research relevant to 
IFE has also been performed, largely for its potential de- 
fense applications. As a result, many MFE and IFE con- 
cepts generally enjoy a far more advanced knowledge 
base than other fusion concepts such as the colliding 
beam and inertial electrostatic concepts. Past efforts 
have been much less extensive both in theory and ex- 
periment, and knowledge about the prospects is far 
more speculative.The likelihood that some alternate con- 
cept may attain and exceed the expected technical and 
economic performance of the tokamak remains specula- 
tive. Developing comparative information judging the 
relative strengths and weakness of a broad range of alter- 
nate concepts and assessing the information base has 
not been a priority of the fusion energy program. In par- 
ticular, there is no current, published DOE-sponsored 
analysis of the comparative technical prospects and chal- 
lenges of the broad array of fusion concepts including 
novel ones or those previously examined and no longer 
pursued. DOE has sponsored and published, however, 
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reviews of alternate MFE concepts that discuss their 
relative level of development and likely prospects, (see 
footnote 10) and has supported some analyses of the 
relative prospects of IFE. (see footnote 11) The lack of 
comparative assessment of non-MFE or IFE concepts is 
consistent with the fusion energy program's primary fo- 
cus on MFE concepts rather than a broader array of fu- 
sion concepts. 



Some alternate MFE concepts previously investigated 
and found less promising than the tokamak may warrant 
reconsideration, based on improvements in technology 
and theoretical understanding. For example, one of the 
major challenges with the stellarator concept was de- 
signing and fabricating the relatively intricate magnets 
required. However, advanced computer-based analytical 
capabilities continue to improve the ability to design 
and manufacture magnets. Some of these techniques 
were developed and used in producing the now prema- 
turely retired Advanced Toroidal Facility (ATF), the 
most recent stellarator. (see footnotel4) While the stella- 
rator may not ultimately prove more attractive than the 
tokamak, improving magnet technology continues to 
reduce one of its principal drawbacks. Advantages rela- 
tive to the tokamak include that they are inherently 
steady state, have no plasma current, and thus do not 
suffer from disruptions and instabilities of the plasma. 
The approximately $1 -billion Large Helical Device 
(LHD), under construction in Japan, is a supercon- 
ducting stellarator similar to ATF in concept, but closer 
to TPX in scope and cost. A similar scale stellarator has 
been proposed in Germany. A much smaller stellarator 
with a cost of about $3 million is under construction at 
the University of Wisconsin as part of DOE's small 
program for alternate fusion concepts. 
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Edge Transport and Turbulence 
Studies on U-3M 

Joint measurements of equilibrium plasma potential, 
electron temperature and ion density, and the fluctua- 
tions in plasma potential and ion density are being car- 
ried out in the edge region of the I = 3 m = 9 
Uragan-3M (U-3M) torsatron. In particular, this work is 
undertaken to clarify the effect of low-frequency electro- 
static turbulence on particle and energy transport. 

In view of data summarized in Ref. 1, the elucidation of 
this question under U-3M conditions is of special inter- 
est because in U-3M the plasma is heated by RF fields 
in the multimode Alfven resonance regime (CO < C0 ci ) 
[2]. It follows from both the theory and some experi- 
mental work that under these conditions parametric in- 
teraction can arise in the edge plasma between the pump 
oscillations and natural electrostatic ion cyclotron har- 
monic oscillations; this can result in an increase of low- 
frequency (drift) turbulence. In this light it would be of 




Fig. 1 . Turbulence-driven particle flow as a function of 
minor radius (the boundary is at ~ 12.5 cm). 
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interest to estimate the equilibrium particle and heat 
flows at the boundary of the confinement region and to 
compare them with the turbulence-induced flows near 
and far from the RF antenna. 

At present the measurements have been made far from 
the antenna in so-called "low-density regime" 
(f RF = 200 kW, B = 0.45 T, n £ « 4X10* 2 cm" 3 , 
r e (0) = 300 eV, r.(0) = 150 eV) [2]. An array of four 
Langmuir probes is used which is introduced from the 
low field side and moved along the major radius. In the 
absence of a material limiter the measurements of float- 
ing potential carried out near the calculated last closed 
magnetic surface (LCMS, a =12.5 cm) indicate the ra- 
dial electric field E r shear layer, which is ~ 2 cm wide. 
In this region the radial profiles of electron temperature 
and ion saturation current are measured. Supposing that 
there are no charged particle sources in the boundary 
layer far from the antenna and that the balance takes 
place between the perpendicular and parallel particle 
flows in the scrape-off layer (SOL), the estimation of 
the equilibrium perpendicular diffusion flow through 
the boundary of the confinement region yields 
r ~ 6 x lO 1 ^ cm 2 s l . This value is 2 to 3 orders of 
magnitude larger than that following from neoclassical 
predictions, but is comparable with T eq values having 
been obtained by similar technique in other closed mag- 
netic traps [1] (see T SOL in Fig. 2). 

At the boundary of the confinement region the relative 
levels of plasma potential (floating potential) and ion 
density (ion saturation current) fluctuations amount, re- 
spectively, to eV /T < 1 and n/ n~ = 0.3. 

Similar to the Texas Experimental Tokamak and the 
ATF torsatron cases, a mean phase velocity shear layer 
exists in the vicinity of the last closed magnetic surface. 
The maximum radial particle flow induced by the turbu- 
lence attains T ~ 5 x 10 16 cm~ 2 s _1 as is shown in Fig. 1 . 
The x coordinate is the distance from the center of the 
poloidal cross-section which is close to the local equilib- 
rium flow. Note that there exists a radial region with T 
going inward outside the confinement region. This 
might be associated with ionization-driven drift wave 
turbulence [3] ; though to confirm or to refute this state- 
ment further, more detailed studies are necessary. 

In Fig. 2, the result of a comparison of the equilibrium 
and turbulence-induced particle flows is presented to- 
gether with the data on similar comparisons made on 
some other closed magnetic traps [1]. We may conclude 
that similar to some tokamaks and the ATF torsatron, in 
U-3M the perpendicular particle flow at the edge is 
dominated by the turbulence-induced flow. In their abso- 
lute values, at least far from the RF antenna, these flows 
are comparable with the edge flows in other closed mag- 



netic traps. In this sense we may say that the chosen 
method of plasma heating in U-3M (multimode Alfven 
resonance) does not result in any significant deteriora- 
tion of particle confinement in comparison with other 
methods of plasma heating (ohmic, electron cyclotron 
resonance, neutral beam injection). 

Victor V. Chechkin for U-3M group 

Institute of Plasma Physics 

National Science Center 

Kharkov Institute of Physics and Technology 

Kharkov, Ukraine 

FAX:007-057-235-2664 
E-mail:ipp@ipp.kharkov.ua 
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Fig. 2. A comparison between turbulence-driven particle 
flow and equilibrium radial flow from Ref. [1] with the 
U-3M data added. 
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Topological aspects of island 
divertor studies on W7-AS 

W7-AS exhibits a flexible boundary topology that can 
be controlled by means of the external rotational trans- 
form. The structure of the plasma boundary is charac- 
terized by the presence of 51m natural islands, and this 
topology can be used to study the feasibility of island 
divertors for power and particle exhaust. The 51m is- 
lands can exist as a closed chain, or they can be inter- 
sected by the inboard limiters. These limiters provide a 
nearly homogeneous scrape -off layer and can be used 
for preliminary investigations to explore the divertor 
potential of W7-AS. This is especially important in con- 
text of the planned island divertor for W7-X. 

The viability of the island divertor concept depends criti- 
cally on the stability of the island structure against exter- 
nal and internal magnetic field perturbations. To assess 
the relevance of vacuum field structures and of perturba- 
tive fields for plasma transport and recycling at the 
edge, we have investigated the boundary topology for 
different values of the edge rotational transform+ a 
within the range of* between 0.34 and 0.64. The ex- 
periments were carried out for net current-free, low-beta 
(~ 0.1%), flat-top electron cyclotron resonance heating 
(ECRH) discharges with P heat = 168 kW and <n £ > = 
3.8 X 10 18 m~ 3 . A Langmuir probe array located at a 



separatrix 
symmetry: 
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Fig. 1 . Tangential view of the plasma taken with a video 
camera installed in Module 5. 



fixed radius outside the separatrix and covering about 
one-third of the poloidal plasma circumference meas- 
ured the poloidal « e ( and T e ) profiles at the boundary, 
while the distribution of the neutral gas density was 
monitored by measuring the spatial variation of the H a 
intensity. Video cameras were used as monitors for the 
visible radiation at the edge. Experimental results from 
these diagnostics show clear evidence of the presence of 
the 5/8, 5/9, 5/10 and 5/11 island chains at the plasma 
edge. Figure 1 shows a tangential view of the plasma 
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TILE INDEX 

Fig. 3. Poloidal profile of the H a intensity in front of the lower inboard limiter. The vacuum field configuration is also 
shown. 



and the upper inboard limiter in Module 5 for a dis- 
charge with+ a = 0.525. The 5/10 island structure at the 
boundary is clearly visible in the picture. 

Figure 2 shows n g profiles at the edge as a function of 
+ . The higher density, higher temperature areas around 
the closed islands are indicated by the bright areas in the 
figure. The change in the separatrix symmetry between 
odd and even values of m (X-point at the midplane for 
5/8 and 5/10 islands, O-point for 5/9 and 5/1 1 islands) 
can be clearly seen in the n g profiles. As+ a decreases the 
islands shift radially outward and are eventually inter- 
sected by the inboard limiters. This leads to a loss of 
confinement within the islands, resulting in the decay of 
the plasma pressure within the island core. The sharp 
edges of the bright areas in Fig. 2 arise due to this sud- 
den loss in confinement. The transition from closed to 
open islands is reflected in the experimental data as a 
doubling of the maxima in the temperature and density 
profiles (best resolved for+ a = 5/9 and 5/10). This is be- 
cause the intersection with the limiters results in the for- 
mation of two SOL branches, and the interaction with 
the limiter now occurs along two zones per island. The 
figure also demonstrates that the 51m structures are not 
significantly perturbed by lower order resonances even 
at+ a = 1/2, where such perturbations are expected to 
have a major influence. Measurements carried out at 
higher beta ((3 = 0.8%) indicate that despite deviations 
in detail, the 51m island chains are still dominant under 
these conditions. 

Comparison of the experimentally measured poloidal 
phase of the islands, in particular their location with re- 
spect to the inboard limiters, with that predicted by vac- 
uum field calculations is a good measure of the 
relevance of these structures for a real plasma. Figure 3 



shows a poloidal H a intensity profile along the surface 
of the lower inboard limiter in Module 3 along with the 
corresponding magnetic configuration. The edge rota- 
tional transform was fixed at 0.510 for this discharge. It 
is clearly seen that the predicted positions of the islands 
(at limiter tiles 9, 13 and 16) are in excellent agreement 
with the location of the maxima in the H a profile. This 
is true also for the 5/8, 5/9, and 5/11 island chains. Poloi- 
dally resolved calorimetric measurements on the in- 
board limiters also confirm these results. 

J. Das for the W7-ASTeam 

Max Planck Institut fur Plasmaphysik 

IPP Garching, Germany 

Phone:0049-89-3299-1725 
FAX: 0049-89-3299-2584 
E-mail: jad@ipp-garching.mpg.de 
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140-GHZ second harmonic 
O-mode electron cyclotron 
heating in W7-AS 

Electron cyclotron resonance heating (ECRH) at the sec- 
ond harmonic ordinary mode polarization (02-mode) 
has a low single -pass absorption at the typical W7-AS 
plasma parameters but gives access to a density twice 
the cut-off density limit of the second harmonic extraor- 
dinary mode (X2-mode). For W7-AS therefore, X2- 
mode heating with its high single-pass absorption is the 
standard heating scenario. In the next step stellarator 
W7-X, a much higher electron temperature and a much 
higher single-pass absorption of the 02-mode are ex- 
pected; therefore the 02-mode ECRH is an attractive 
scenario [1]. With this in mind we have investigated the 
general absorption properties of the 02-mode in W7-AS. 

A 140-GHz, 40-ms pulse with 0.7 MW power in 02- 
polarization was launched at an optimum oblique angle 
of 20° into a NBI-sustained target plasma with a mag- 
netic field of 2.5 T. A central electron density of 1 .8 x 

20 —3 

10 m was chosen, which is well above the X2 cut- 
off density. The electron temperature is only 0.4 keV un- 
der these high-density conditions. The plasma energy 
content represented by the diamagnetic signal in Fig. 1 
increased by about 1 .5 kJ compared to a similar dis- 
charge with NBI only. The absorbed heating power is 
110 kW, determined from the time derivative of the dia- 
magnetic signal at the heating pulse switch-off. This is 
roughly twice the calculated power absorbed in a single 
pass and can be explained by wave reflection at the in- 
ner torus wall after a first plasma transit. 

The soft X-ray emission was analyzed to determine the 
power deposition zone. The line-integrated X-ray emis- 
sions along nine different chords more or less in the di- 
rection of the plasma elongation are recorded with a 
spatial resolution of about 2 cm. Despite the low spatial 
resolution the time delay between the decay of the dif- 
ferent X-ray detector signals after heating power switch- 
off is easily resolved. This time delay with respect to the 
moment of power switch-off is sketched in Fig. 2 as a 
function of the radial position of the X-ray emission. 
The time delay increases with the effective chord radius, 
indicating that the heating power must have been ab- 
sorbed in the dense plasma core at the ECR layer. 

H. Laqua for the W7-AS Team 

Max-Planck-lnstitut fur Plasmaphysik, EURATOM Ass. 

D-85748 Garching, FRG 

E-mail laqua@ipp-garching.mpg.de 
Phone:0049-89-3299-1963 
FAX: 0049-89-3299-2584 
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Fig. 1 . Energy content from the diamagnetic signal of a 
discharge with combined NBI and 02-mode ECRH (up- 
per curve, suppressed origin). For comparison a purely 
NBI-heated discharge is shown in the lower curve. The 
time window of the ECRH pulse is shown by the mark- 
ers. 
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Equilibrium plasma currents in 
quasi-symmetric stellarators with 
a helical magnetic axis 

The concept of a quasi-symmetric stellarator was ad- 
vanced in [1,2]. A latent symmetry of the three- 
dimensional (3-D) magnetic field in such systems must 
improve their confinement properties because of the par- 
ticle drifts and heat fluxes across the magnetic surfaces. 
The finite-beta plasma currents can also be reduced. In- 
itially, such configurations were investigated with the 
use of 3-D equilibrium codes, which can require signifi- 
cant computer resources. For example, the Wendelstein 
7-X configuration was found computationally. 

More recently, theoretical and numerical methods have 
been developed to make examination of some of the 
questions related to quasi-symmetry (QS) easier. Theo- 
retical investigations show [3-5] that in the paraxial ap- 
proximation the QS conditions are given as relations of 
following parameters: k, K, e and 8, where k(s) and k(^) 
are the first and the second curvatures of the magnetic 
axis, E(s) is the ellipticity of the magnetic surface cross- 
sections, 8(s) is the angle between the normal to the 
magnetic axis and the small semi axis of the ellipse, and 
s is the arc length of the magnetic axis. The correlation 
equations for these parameters were obtained (in the 
first order of the paraxial expansion) in [3]. 

In this work the equilibrium plasma currents are investi- 
gated in the paraxial approximation for quasi-symmetric 
configurations with a helical magnetic axis where QS 
conditions are determined by the equations obtained in 
[3]. The parameter 



(1) 



which is one of the important stellarator optimization 
parameters, is also calculated. In Eq. (1) j n and j ± are the 
plasma current components, parallel and perpendicular 
to the magnetic field, and the integrals imply integration 
over the volume of the layer between closely spaced 
magnetic surfaces. 

To solve the equations [3] for the QS conditions, nu- 
merical methods are necessary. Because of the condi- 
tions the magnetic surface geometry in the systems 
under consideration is very complicated. Therefore, the 
j'll and j ± and the geometric factor [Eq. (1)] are calcu- 
lated by numerical integration along the field lines using 
the methods of paper [6]. 

We start from an assignment of the magnetic axis shape 
and find the parameters e and 8 as functions of 5 using 



the equations of paper [3]. After that we obtain the ex- 
pressions for the magnetic field, B, which are necessary 
for the plasma current differential equations. For this 
purpose we use the work coordinate system (a, 9, s), 

-ry 

x = ae 



cos6 = p cos(i3+8) , 
= a e^ 2 sin6 = p «'n(i3+8) , 



(2) 



where E(s) = tanhr|(s); a is the magnetic surface label, 9 
is a quasi-poloidal angle, p is the distance from the mag- 
netic axis, and •& is the angle measured from the normal 
to the magnetic axis. The expression for B may be repre- 
sented as in [3] (on the magnetic axis, the magnetic 
field, B Q , is constant, and we assume that the total 
plasma current, 7, is zero): 



2nB = FV^ + V^, 



or 



B 



F_ 
2% 



(3) 



(4) 



Here V = V - Va(VaV)/IVal, F=2%RB Q , R = L/2%, L 
is the full length of the magnetic axis, £ = s/ R, a/ R is 
the small parameter of the paraxial expansion, 

1 3(0 a 2 r (dd \ „„ coshri dr\ . „ -i 
1 ! A — -° cos29 + , ' ... sin29 



f39 




2 dt; 



(5) 



In the approximation considered, the term (3(p/3Q/F 
may be omitted in comparison with unity. The Va, V9 
and V£ components can be calculated using the expres- 
sions in Eq. (2). 

In accordance with the calculation methods developed 
in [6] the currents j., and j ± and the parameter [Eq. (1)] 
can be calculated using the system of ordinary differen- 
tial equations. It includes first, the equations of the mag- 
netic field lines 



da/dx = 0, dQ/dx = BVQ, ds/dx = BVs. 

Second the equations for the equilibrium plasma cur- 
rents are given as 

j L = cp'[BVY\/B 2 j n = cp'hB, 



dh/dx = -2[BV5]V V P/5 J 



(6) 



(7) 



(8) 



where c is the speed of light, and h arises from the equi- 
librium equations. Third, the equations for calculating 
the integrals of these currents over the volume of the 
layer between neighboring magnetic surfaces follow 

df/dx = (hBB Q ) 2 , 
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df 2 /dx = {B ^/B) 2 f 1 (0)=f 2 (0) = 



0) 



(the parameter x is an integration variable, the prime de- 
notes the derivative with respect to the magnetic surface 
function «P (BV*P = 0), the definition V*P = Va/IVa l is 
used, Va Q is the Va value at the starting point of the in- 
tegration). The right-hand sides of Eq. (9) are propor- 
tional to the squares of j n and j ± . Therefore, the 
parameter [Eq. (1)] is determined as 



X — . 

m 



(10) 



X m is the upper limit of the integration interval. 

We considered a magnetic axis specified by the paramet- 
ric form 



r = r Q + rjCos(M])), z = rjSin(M|)), 



(11) 



which is a helical line on a torus (r, (]), z are the cylindri- 
cal coordinates, N is a number of periods of the system). 
The calculations were carried out for several configura- 
tions at r/r = 0. 132, N = 5 and N = 3. The results of 
the calculations of the factor Eq. (1) (for a/r l = 1/13.2 
at N = 5 and a/r l = 1/26.4 at N = 3) are presented in 
Table 1, where the parameters E Q , n,, \i, T|, and 8, which 
characterize the configurations under consideration, are 
also presented. Here E = e* 1 ^ is the ratio of large and 
small semiaxes of ellipse at (j) = 0, n, = nIN; n is the 
number of turns made by the ellipse relative to the nor- 
mal to the magnetic axis on the full length of the axis; (J, 
is the rotational transform relative to the normal to the 
magnetic axis; T| j 2 3 and 8j 2 3 mean the corresponding 
distributions of T) and 8 over s, which are determined by 
the assignment of E Q and 8(0) = 0. The parameters in 
the (TJj, 8j) and (T| 2 , 8 2 ) cases correspond to the mag- 
netic configurations considered in Ref. [5]. 

Some of the results pertain when the QS conditions are 
violated. The 8 = 8,, T|= T| 2 case corresponds to slightly 
violated conditions, and the 8 = 0, T) = cases corre- 
spond to strongly violated conditions and to a circular 
shape of the magnetic surface cross sections. 

It follows from the results obtained that for the consid- 
ered variants of the systems, the geometric factor 
[Eq. (1)] is less than unity, satisfying one of the impor- 
tant stellarator optimization conditions. For the same 
magnetic axis parameters but at strong violation of the 
QS conditions (for the circular shape of the magnetic 
surface cross sections) the mentioned factor is essen- 
tially greater than unity. 

The decrease of the parameter [Eq. (1)] corresponds to a 
lowering of the secondary plasma currents and to a re- 
duction of distortion of the magnetic configuration due 



to finite plasma pressure. Moreover, in the Pfirsch- 
Schliiter regime, the geometric factor of the diffusion co- 
efficient increase due to the particle drift in the 
inhomogeneous magnetic field is determined by the pa- 
rameter [Eq. (1)]. Thus, using relatively simple methods 
of investigation, we have demonstrated the possibility 
of improvements in these parameters for quasi-symmet- 
ric systems in the Pfirsch-Schliiter regime. 

The results of calculations of the plasma currents were 
also used for the analysis of the QS condition in the in- 
variant form [5] for the systems under consideration. 
The analysis has shown that the QS condition in the in- 
variant form is well satisfied at sufficiently small dis- 
tances from the magnetic axis for the configuration 
variants, for which the QS conditions were obtained in 
the paraxial approximation with the help of the equa- 
tions of paper [3] . 

The work was carried out in the Institute of Plasma 
Physics, National Science Center "Institute of Physics 
and Technology," Kharkov, in collaboration with Rus- 
sian Science Center "Kurchatov Institute." 
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